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ABSTRACT 

Context. Recent VLBI studies of the morphology and kinematics of individual BL Lac objects (S5 1803+784, PKS 0735+178, etc.) 
have revealed a new paradigm for the pc-scale jet kinematics of these sources. Unlike the apparent superluminal outward motions 
usually observed in blazars, most, if not all, jet components in these sources appear to be stationary with respect to the core, while 
exhibiting strong changes in their position angles. As a result, the jet ridge lines of these sources evolve substantially, at times forming 
a wide channel-flow. 

Aims. We investigate the Caltech-Jodrell Bank flat-spectrum (CJF) sample of radio-loud active galaxies to study this new kinematic 
scenario for flat-spectrum AGN. Comparing BL Lac objects and quasars in the CJF, we look for differences with respect to the 
kinematics and morphology of their jet ridge lines. The large number of sources in the CJF sample, together with the excellent 
kinematic data available, allow us to perform a robust statistical analysis in that context. 

Methods. We develop a number of tools that extract information about the apparent linear and angular evolution of the CJF jet ridge 
lines, as well as their morphology. In this way, we study both radial and non-radial apparent motions in the CJF jets. A statistical 
analysis of the extracted information allows us to test this new kinematic scenario and assess the relative importance of non-radial 
and radial motions in flat-spectrum AGN jets, compared to those of quasars. We also use these tools to check the kinematics for 
(multi-wavelength) variable AGN. 

Results. We find that approximately half of the sample shows appreciable apparent jet widths (> lOdegrees), with BL Lac jet ridge 
lines showing significantly larger apparent widths than both quasars and radio galaxies. In addition, BL Lac jet ridge lines are found 
to change their apparent width more strongly. Finally, BL Lac jet ridge lines show the least apparent linear evolution, which translates 
to the smallest apparent expansion speeds for their components. We find compelling evidence supporting a substantially different 
kinematic scenario for flat-spectrum radio-AGN jets and in particular for BL Lac objects. In addition, we find that variability is 
closely related to the properties of a source's jet ridge line. Variable quasars are found to show "BL Lac like" behavior, compared to 
their non-variable counterparts. 

Key words. Galaxies: statistics - Galaxies: active - Galaxies: nuclei - Galaxies: jets - BL Lacertae objects: general 



1. Introduction 

Although observed in the minority of active galaxies (~ 5-15%; 
Kellermann et al. 1989; Padovani 1993; Jiang et al. 2007), ex- 
tragalactic jets are some of the most pronounced moiphological 
features in AGN research. Their presumably direct connection 
to the active core and the supermassive black hole (SMBH) re- 
siding there, makes them invaluable tools in the effort to char- 
acterize the properties and the underlying physics of activity in 
galaxies. VLBI observations enable the direct imaging of AGN 
jets and thus the study of their properties on parsec scales. One 
of the most prominent discoveries, related to jet kinematics, was 
that of apparent superluminal motion of jet components (e.g., 
Whitney et al. 1971; Pearson & Readhead 1981), a combina- 
tion of relativistic expansion speeds (close to the speed of Ught) 
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and the projected geometry onto the plane of the sky (e.g., Rees 
1966). 

Jet kinematics, as studied through the investigation of dis- 
tinct components, is usually explained in terms of the shock-in- 
jet model (e.g., Marscher & Gear 1985), where the observed jet 
knots are manifestations of shocks propagating down the jet at 
relativistic speeds. Beaming and projection effects regulate the 
observed properties of the jets. There has been continuous effort 
to distinguish whether the different types of active galaxies (e.g., 
quasars, BL Lacs, FanaiTof-Riley, etc.) are just a result of ori- 
entation effects, or if additionally these objects have intrinsically 
different properties. The current paradigm is that indeed different 
jet properties can be attributed to geometrical effects combined 
with factors such as the black hole mass or the accretion rate. For 
example, Ghisellini et al. (1993) for a sample of 39 superluminal 
sources find no appreciable difference between the distiibution 
of Doppler factors between BL Lacs and flat-spectrum radio- 
quasars (FSRQs), with RGs showing smaller values. Some indi- 
cations to the contrary also exist (e.g., Gabuzda 1995; Gabuzda 
et al. 2000). Analysis of statistically important samples (large 
number of sources and/or stringent selection criteria) of active 
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galaxies have been of fundamental importance to this end (e.g., 
Ghisellini 1993; Vermeulen & Cohen 1994; Vermeulen 1995; 
Taylor et al. 1996; Hough et al. 2002; Lister & Homan 2005; 
Britzen et al. 2007b). 

Although it is crucial to pursue a statistical approach to the 
open problem of AGN jet kinematics, the study of individual 
sources is indispensable. It has helped to elucidate particular 
mechanisms or effects that might get smoothed out by the poor 
temporal or spatial resolution usually available in big statistical 
samples. Such a case of a detailed study of an individual object 
is that of S5 1803+784. 

1.1. S5 1803+784: A case study 

S5 1803+784 is an active galaxy at a moderate redshift of z=0.68 
(Hewitt & Burbidge 1989). It has been classified as a BL Lac ob- 
ject. Being a member of the complete S5 sample (Witzel 1987) 
it has been extensively studied in the radio, at different wave- 
lengths and with different instruments (see Britzen et al. 2010a 
for a detailed recounting of the source's radio observations his- 
tory). Typical for this class of objects, 1803+784 has been ob- 
served to be variable in the radio and the optical on both long and 
short-timescales (e.g., Wagner & Witzel 1995; Heidt & Wagner 
1996; Nesci et al. 2002; Fan et al. 2007). 

Britzen et al. (2010a), by (re-)analyzing more than 90 epochs 
of global VLBI and VLBA data, reveal a new kinematic scheme 
for 1803+784. All components in the inner part of the jet (up 
to 12 mas) appear to remain stationary with respect to the core. 
This behavior is seen at all frequencies studied by the authors 
(1.6 - 15 GHz). In contrast to this, the components show strong 
changes in their position angles, implying a prevailing move- 
ment perpendicular to the jet axis. 

Britzen et al. (2010a) also studied the jet ridge line morphol- 
ogy and evolution of S5 1803+784. A jet ridge line at a given 
epoch is defined as the line that linearly connects all compo- 
nent positions at that epoch. The authors find that the jet ridge 
line changes in an almost periodic manner, starting resembling a 
straight line, evolving into a sinusoid-like pattern, and finally re- 
turning to its original linear pattern, although slightly displaced 
from its original position. A period of ~ 8.5 years is calculated 
for the evolution of the jet ridge line. 

Finally, the authors find that the jet changes its apparent 
width (in a range between a few and a few tens of degrees) in 
an almost periodic way with a timescale similar to the one found 
from the evolution of the jet ridge line. All of the above proper- 
ties support a new kinematic scheme for 1803+784, where com- 
ponents follow oscillatory-like trajectories, with their movement 
predominantly happening perpendicular to the jet axis rather 
than along it. Moreover, the jet appears at times to form a wide 
channel of flow, while changing its width considerably across 
time. 

1.2. Motivation 

S5 1803+784 is one of several BL objects exhibiting the behav- 
ior described above. 0716+714 has been shown to behave in a 
similar way (Britzen et al. 2009), with most of its components 
being stationary with respect to the core while changing their 
position angles considerably. PKS 0735+178 is another exam- 
ple of a source with similar, but rather more complicated, kine- 
matic properties (Gomez et al. 2001; Agudo et al. 2006; Britzen 
et al. 2010b). Under the unification scheme of active galaxies 
(e.g., Antonucci 1993; Urry & Padovani 1995), both BL Lac ob- 



jects and FSRQs are believed to be active galaxies for which the 
viewing angle to their jets is very small, leading to strong rela- 
tivistic effects. Given the similar viewing angle distributions, in 
the last several years the general classification of blazar has often 
been used to describe members of either class, also in terms of 
their jet properties. However, this phenomenological unification 
of FSRQs and BL Lacs comes into question in light of the recent 
investigation of sources like S5 1803+784, PKS 0735+178 and 
0716+714. Is the peculiar kinematic behavior seen in these ob- 
jects revealed due to the unprecedented richness of the datasets 
available, and could therefore be relevant for all flat-spectrum 
radio- AGN, or are BL Lacs characterized by a genuinely differ- 
ent set of kinematic properties? 

This paper focuses on the relevance of this new kinematic 
scheme for flat-spectrum radio-AGN, while investigating the ap- 
parent divide between BL Lac and FSRQ jet kinematics. As 
valuable as single source studies are to an in-depth understand- 
ing of particular phenomena, there are a number of biases or 
unaccounted factors that alter and ultimately hinder a universal 
application of their results. In this context, we use the CJF sam- 
ple to statistically investigate and assess the similarity, or diver- 
gence, of the kinematic and morphological properties between 
the two distinct sub-samples of FSRQs and BL Lac objects in the 
CJF. We want to test whether jet components of BL Lac objects 
indeed show slower apparent speeds with respect to their cores 
compared to FSRQs. Furthermore we are interested in the phe- 
nomenon exhibited in S5 1803+784 of an, at times, very wide 
jet, as well as a strong evolution of that width. For this investi- 
gation we use tools that extract information from the jet ridge 
line of the sources, instead of focusing on individual compo- 
nents. This allows for a investigation that is mostly independent 
of component modeling and cross-identification. 

The paper is organized as follows: in Sect. 1.1 we introduce 
the new kinematic scheme for BL Lac objects, as shown for the 
case of the source 1803+784, and discuss the motivation of this 
work, in Sect. 2 we describe the CJF sample, in Sect. 3 we de- 
scribe the data used, in Sect. 4 we present the analysis of our data 
and the results, and in Sect. 5 we discuss our results and give 
some conclusions. Throughout the paper, we assume the cos- 
mological parameters Hq = 71 kms^^Mpc^^, Q.m - 0.27, and 
Qa = 0.73 (from the first-year WMAP observations; Spergel 
et al. 2003). 



2. The CJF sample 

The CJF sample (Taylor et al. 1996) consists of 293 radio-loud 
active galaxies selected (see Table 1) from three different sam- 
ples (for details, see Britzen et al. 2007b). The sources span a 
large redshift range (see Fig. 1), the farthest object being at a 
redshift z = 3.889 (1745+624; Hook et al. 1995) and the clos- 
est at z = 0.0108 (1146+596; de Vaucouleurs 1991). The aver- 
age redshift of the sample is Zavg - 1.254, ZsLUicavg - 0.546, 
ZRCavg - 0.554, and ZFSRQ,m-g - 1.489 for BL Lacs, radio 
galaxies, and FSRQs, respectively. All the objects have been ob- 
served with the VLBA and/or the global VLBI network. Each 
source has at least 3 epochs of observations (with a maximum 
of 5 epochs) and has been imaged and studied kinematically 
(Britzen et al. 1999; Britzen et al. 2007b; Britzen et al. 2008). 
The X-ray properties have been studied and correlated with their 
VLBI properties (Britzen et al. 2007a). The evolution of active 
galaxies, in the context of the merger-driven evolution scheme 
(e.g., Hopkins et al. 2006), has also been investigated with the 
help of the CJF, identifying candidate CJF sources in different 
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Fig. 1. Redshift distribution of radio galaxies (grey blocks), BL 
Lacs (dashed line), and FSRQs (solid line) in the CJF sample. 



evolutionary stages, including new binary black hole candidates 
(Karouzos et al. 2010). 



Table 1. CJF sample and its properties. 



Frequency(MHz) 


4850 


Flux lower limit @5GHz 


350mJy 


Spectral Index 


4850 > _Q 


Declination 


5> 35° 


Galactic latitude 


\b\ > 10° 


# Quasars 


198 


# BL Lac 


32 


# Radio Galaxies 


52 


# Unclassified 


11 


# Total 


293 



2.1. Radio emission 

The CJF is a flux-limited radio-selected sample of flat-spectrum 
radio-loud AGN. The sample was originally created to study, 
among other things, the kinematics of pc-scale jets and apparent 
superluminal motion (see Taylor et al. 1996 for more details). 

The CJF sample (Table 1) has been most extensively stud- 
ied in the radio regime (e.g., Taylor et al. 1996; Pearson et al. 
1998; Britzen et al. 1999; Vermeulen et al. 2003; Pollack et al. 
2003; Lowe et al. 2007; Britzen et al. 2007b; Britzen et al. 2008). 
Britzen et al. (2008) developed a localized method for calculat- 
ing the bending of the jet associated with individual components. 
The maximum of the distribution of local angles is at zero de- 
grees, although a substantial fraction shows some bending (0-40 
degrees). A few sources exhibit sharp bends of the order of > 50 
degrees (see Fig. 13 in Britzen et al. 2008). 

Although the CJF sample consists mostly of core-dominated 
AGN, presumably highly beamed sources, the kinematical study 
of the sample identifies a large number of sources with station- 
ary, subluminal, or, at best, mildly superluminal outward veloc- 
ities (e.g., see Fig. 15 in Britzen et al. 2008). Combined with 
a number of sources with inwardly moving components (e.g., 
0600H-422, 1751-H441, 1543-H517, Britzen et al. 2007b), these 
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sources do not fit into the regular paradigm of outward, superlu- 
minaly moving components in blazar jets. One explanation of 
these peculiar kinematic behaviors is that of a precessing, or 
helical jet (e.g., Conway & Murphy 1993) possibly as a result 
of a SMBH binary system. Other interpretations include trailing 
shocks produced in the wake of a single perturbation propagat- 
ing down the jet (e.g., Agudo et al. 2001; but also see Mimica 
et al. 2009), or standing re-collimation shocks (e.g., Gomez et al. 
1995). 

Karouzos et al. (2010) compile a list of all the CJF sources 
that have been found to show long timescale variability in the ra- 
dio (as well as in other wavelength regimes). This list comprises 
in total 40 CJF sources, 27 of which have been argued to show 
possibly periodical variability of their fluxes. The authors do not 
take into account intra-day variability. 

3. Data 

The work presented here is heavily based on the kinematic anal- 
ysis of the CJF sample (Britzen et al. 2007b; Britzen et al. 2008). 
An extensive observing campaign of all 293 CJF sources was un- 
dertaken using both the VLBA and the global VLBI array at 5 
GHz (see Britzen et al. 2007b for details). We note that five CJF 
sources were initially excluded from any further analysis due 
to problematic observations (0256-h424, 0344-H405, 0424-H670, 
0945+664, 1545-H497) and therefore are also omitted here. In 
total, 288 sources are considered and analyzed in the following 
sections. Of these, according to the optical classification from 
Britzen et al. (2007b), 196 are classified as quasars, 49 as radio 
galaxies, 33 as BL Lac objects , and 10 are not classified. 

Due to the scope of the CJF program, the identification and 
analysis of pc-scale jet component kinematics has focused on 
the part of the jet that is beamed towards us. For a number of 
sources, several components belonging to the counter-jet have 
been identified. However, for these sources cross-identification 
of the counter-jet components over epochs has not been carried 
out. For this reason, and given the nature of the analysis that we 
undertook (see below), we have excluded all counter-jet compo- 
nents in the following investigation. In a total number of 2468 
components identified, 82 (3.32%) counter-jet components have 
been identified and excluded from our analysis. Britzen et al. 
(2007b) report that on average radio galaxies have 3.6 compo- 
nents identified per jet, 2.7 components are identified per quasar 
jet, and 2.9 components per BL Lac jet. This reflects the relative 
difference of projected jet length for the different types of object 
classes. This shall be discussed more thoroughly in Sect. 5. 

In the following, the tools that we use for the analysis of 
the CJF jet ridge lines, in the context described in Sect. 1.1, are 
described. In short we used the following measures: 

- Monotonicity Index, M.I. 

- Apparent Jet Width, dP 

- Apparent Jet Width Evolution, AP 

- Apparent Jet Linear Evolution, A{ 

We note that all the above measures, as their names imply, refer 
strictly to values projected onto the plane of the sky. Although it 
is possible to constrain, or in some cases have a specific estimate 
of the viewing angle of each source and therefore attempt to cal- 
culate intrinsic jet properties, this is outside the scope of this pa- 
per For the following sections, we adopt the basic assumption as 
described by the AGN unification scheme (e.g., Antonucci 1993, 
Urry & Padovani 1995) that BL Lacs and FSRQs are seen at the 
smallest viewing angles, while radio galaxies have their jet axis 
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further away from our line of sight. As we are mainly interested 
in the comparison between FSRQs and BL Lacs, deprojection 
of the jet properties investigated here is not critical. For the fol- 
lowing analysis, and for the sake of brevity, we shall drop the 
characterization of apparent for each of these values, although 
this will be implied throughout unless otherwise stated. 

3.1. Monotonicity Index, M.I. 

It is known that a large number of active galaxies exhibit bent 
or otherwise non-linear jet morphologies on various scales. 
Individual sources like S5 1803+784 and PKS 0735+178 (see 
references above), as well as others (e.g. 3C 345; Lobanov & 
Zensus 1999; B0605-085; Kudryavtseva et al. 2010), have been 
extensively studied to understand the origin of these bends. 
Britzen et al. (2008) calculate the "bends" between successive 
components of the pc-scale jet for all CJF sources, finding large 
local variations of the position angles of the sources. In this con- 
text, we are interested in quantifying the bending of the whole 
jet ridge line. Moreover, we want to differentiate between a 
"monotonically-bent" jet, i.e., a jet that is bent in only one di- 
rection (see Fig. 2, right), as opposed to a more sinusoid-like 
morphology (similar to what is seen for some epochs of S5 
1803+784, see Fig. 2, left). This is done by means of the mono- 
tonicity index, M.I.. 

We quantify a sinusoid-like morphology of a jet by identify- 
ing the local extrema in a given jet ridge line, in the core sep- 
aration - position angle plane. For an epoch /, a component m 
exhibits a local extremum under the definition 

\0m - ()in±l\ ^ I0(d6,„ + d6,„±i), 

where 0,„ and d9„, denote the position angle of component and 
its uncertainty. 0,„ is calculated as 

X 

6„, - arctan — , 

where X and Y are cartesian coordinates on the plane of the sky. 
Having calculated the number of extrema for a given jet ridge 
line at an epoch i, we define the M.I. as 

number of extrema 



where N is the total number of components at that given epoch. 
This is a crude calculation, but can give us a handle on how the 
bending of the jet behaves along the jet. For M.I. values close to 
one, the jet resembles more a sinusoid. An M.I. value close to 
zero reveals a monotonic, single-bend, jet morphology. We nor- 
malize for the number of components N to account for longer, or 
shorter, jets and to enable comparison between different sources. 
As an example, for the sources shown in Fig. 2, 1928+738 is 
found to have an M.I. value of 0.4, compared to 0700+470, that 
gives an M.I. value of 0. 

As is the case for most of the tools described in this sec- 
tion, the value of M.I. depends heavily on the resolution of the 
observations. Given that the resolution is not the same across 
all epochs and sources, an uncertainty is introduced when com- 
paring two M.I. values of two different sources. A final remark 
pertains to the definition of an extremum. We use a lOcr' value 

' Here a is defined as the sum of the position angle errors for 6„, 
and 9m±\ ■ The choice of this cr reflects an original underestimation (in 
Britzen et al. 2007b) of the position angle errors during the component 
fitting. 



as the lower limit for flagging an extremum. A different level of 
significance (e.g., 5cr) would result in different values of M.I.. 
The choice of the lOcr significance is a conservative approach to 
the identification of extrema in the jet ridge line. In the follow- 
ing, M.I. shall be used as more of a qualitative tool rather than a 
quantitative measure of the actual jet morphology. 

3.2. Apparent jet width, dP 

For an epoch / and a jet consisting of N components charac- 
terized by their core separation and position angle (r,, 6*,), we 
identify the components with the maximum and the minimum 
position angles. The apparent width of the jet dP, measured in 
degrees, is then calculated as (see Fig. 3): 

dPi = 0™^ - fff", 

while the error is calculated by the propagation of errors for- 
mula. There is a number of factors that should be considered at 
this point. The most obvious drawback for the above definition 
is the non-localized nature of this measure. Using two differ- 
ent components, at different core separations, gives us only an 
approximate notion of the width of the flow. Higher resolution 
investigations might resolve the structure of the jet perpendicu- 
lar to the jet axis and provide us with a more real estimate of the 
local width of the jet. Alternatively, one can calculate a localized 
value of the jet width by using the FWHM of the fitted Gaussian 
components (e.g., Pushkarev et al. 2009). In this case however, 
the width is heavily dependent on the beam size and therefore on 
the resolution. Instead, by using the definition above, we seek to 
quantify the opening of the jet flow and identify the effect of a 
channel-like jet, as seen in the case of S5 1803+784. In the con- 
text of the localization problem, we have implemented in our cal- 
culations a constraint in the maximum allowed distance between 
the two components identified with the minimum and maximum 
position angles. The value of this limit depends on the jet length 
and shall be discussed in Sect. 4. 

It is clear that in the case of a narrow but highly bent jet, dP 
would then measure the amplitude of the bending rather than the 
width of the flow itself. In the context of this study such an effect 
does not present a problem as we are, to first order, interested in 
the final result, i.e., an apparently wide jet ridge line opening, 
rather than whether this effect is due to a genuinely wide jet flow 
or rather a high-amplitude bending. 

3.3. Apparent jet width evolution, AP 

The jet width evolution AP, measured in degrees per unit time, 
is measured between two successive epochs ; and (i-1). It is cal- 
culated as follows: 

dPj - dPj-i 

where T,- denotes the time at epoch /, measured in years. We can 
also define the maximum jet width evolution as; 

AP"'"-' = max{APi,j}. 

This value is characteristic of each source and reflects the max- 
imum potential width change of the jet flow for that source. We 
should note here that the jet width evolution does not reflect an 
angular speed (as implied by the units of degrees per unit time). 
Instead the per unit time reflects a normalization for time that 
ensures a comparison between sources that have been observed 
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Fig. 2. Radio maps at 5GHz of 1928+738 (left) and 0700+470 (right), with jet ridge lines superimposed. Two examples of a 
sinusoid-like jet morphology (left) and a single-bent one (right). VLBI maps from Britzen et al. (2007b). 



during different time spans. The currently available data are too 
sparse to allow us a calculation of a real angular speed for how 
the jet width (and orientation) changes with time. 

For both the width dP and the width evolution AP, the sen- 
sitivity and dynamic range of our observations play a regulating 
role. Given that we need at least two components to define the 
width, a dimming of one of these components across epochs can 
result to a false or altered value for both the calculated width and 
width evolution. A more detailed discussion of this can be found 
in Sect. 5. 

3.4. Apparent jet linear evolution, Kt 

We calculate the linear evolution across all available epochs and 
for all components, ultimately acquiring a value reflecting the to- 
tal linear displacement of the whole jet ridge line. We use plane- 
of-the-sky coordinates (X,, F,) to calculate the linear displace- 
ment of component m between epochs / and { i+1 ): 

If = Vax2 + Ay2. 

To calculate the total displacement of the whole jet ridge line 
we then need to sum up over all components and all available 
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Fig. 3. Radio map of 0700+470 at 5GHz, exhibiting the calcu- 
lation of the jet ridge line width. The two arrows denote the 
components at maximum and minimum position angle, while 
the area in between is defined as the opening, or apparent width, 
of the jet ridge line. The VLBI core is found at (0,0) coordinates. 
VLBI map from Britzen et al. (2007b) 



observing epochs: 

1- 1 m 

We need to account for both the different time span of observa- 
tions, as well as the different number of components. Therefore 
we define the jet linear evolution M, measured in parsecs per 
unit time and per component, as follows: 

N(Ti-Ti) NdT' 

where N here is the total number of components used across 
all the epochs. We need to underline a fundamental differ- 
ence between the way the jet linear evolution. At, is calcu- 
lated, compared to the measures described previously. For the 
calculation of the displacement of an individual component £ 
between two consecutive epochs, and consequently of A£, the 
cross-identification of components across epochs is necessary. 
Therefore, for At, unlike the previously discussed jet ridge line 
characteristic values, the actual identification of components is 
important. 

The reason for normalizing this value over the time span of 
observations is the same as for the jet width evolution. In this 
case however we also need to account for the number of com- 
ponents for each individual jet. We therefore divide the total dis- 
placement of the whole jet ridge line by the actual number of 
components used for the calculation (hence the number of com- 
ponents cross-identified across each pair of observation epochs). 

At essentially reflects the apparent speed distribution of all 
cross-identified components of the jet and therefore represents 
a value characteristic for the whole jet, rather than for any 
individual component. By summing up all components and 
epochs we trade temporal and positional resolution for a univer- 
sal treatment of the entire jet. In this way we can test whether 
the kinematics of BL Lac objects is fundamentally different 
than that of FSRQs while averaging out localized properties of 
individual components. A( can be seen as a mean jet component 
speed, with the difference that it is acquired through averaging 
not only over all components, but additionally over all available 
epochs. Although it certainly reflects a measure of the outward 
motion in BL Lac jets, the calculation of Ai is done in such a 
way, that the potential curvature of the components' trajectories 
is taken into account. This separates A{ from a simple linear 
regression fit to the core-separation versus time diagrams 
usually employed to calculate outward velocities, making it 
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sensitive to non-radial motions, that are otherwise missed. 

Lastly, given the large span of redshifts that the CJF sources 
cover, we must account for the different linear scales probed for 
sources at different distances from us. By comparing sources 
at different redshifts, one implicitly studies a different part of 
their jet. To counter this effect, we explicitly define a maximum 
apparent core separation limit, above which components are 
not included in our investigation. There are however additional 
arguments that indicated the necessity of such a limit, both 
technical as well as physical. It is known that AGN jets appear 
most curved closer to the core, with bending effects becoming 
dumped at larger distances (e.g., Krichbaum et al. 1994, Britzen 
et al. 2000). It should also be noted that at larger separation 
from the reference center of a map (the core), the uncertainties 
involved in the identification and fitting of jet components 
increase. By constraining our analysis in the inner-structure of 
the jets, we partly compensate for this effect. In this way we 
make certain that we probe the same linear scales for all CJF 
sources, albeit with different resolution. The latter effect can be 
counteracted by using redshift bins throughout. Additionally, 
through use of redshift bins, we can probe a possible evolution 
of the studied properties with redshift, if any. Given that the 
separation limit used is an apparent one, a level of uncertainty, 
dependent on the scatter of the viewing angle distribution within 
a certain object class, is expected. 

4. Analysis and results 

As the apparent jet length (expressed by the maximum distance 
of the outermost jet component identified in each source, across 
all available epochs) and its morphology (as expressed by the 
M.I.) are used as a basis for the further analysis, we are going to 
briefly discuss these first. 

4.1. Apparent jet length and morphology 

We find BL Lacs to show shorter jets on parsec scales in median 
values, compared to quasars and radio galaxies (17.3 + 2.7 pc 
compared to 27.1 + 2.4 pc and 43 + 7 pc;~ 3cr level), with their 
average values however not being statistically significantly dif- 
ferent (34+5 pc, compared to 3 1 ±4 pc and 50+6 pc, respectively; 
< 3cr)^. Having calculated the average and median apparent jet 
lengths for the CJF objects, we are able to define a maximum 
apparent core separation limit (for the reasons discussed in Sect. 
3). As we are interested in the behavior of BL Lacs, we use them 
as the basis of our decision. For 1803+784, peculiar kinematics 
of the components is observed for the inner-most part of its jet 
(up to 12 mas), with most robust effects out to 6 mas. For this 
source's redshift (z=0.68), this translates roughly to (82 pc) 41 
pc. Similarly to 1803+784, 0735+178 and 0716+714 also show 
the most prominent evolution of their jet ridge line out to ~ 6 
mas (39.6 pc and 33 pc, for their respective redshifts). Given that 
most of the BL Lacs in our sample is at z < 1, the length scale 
implied by the three sources mentioned above (the only objects 
whose jet ridge line has been investigated in detail) ensures that 
the inner-jet of sources at z< 1 is covered. We adopt therefore 
an apparent core separation limit of 40 parsecs. All components 
at separations larger than 40 parsecs with respect to the core are 

^ Thiroughiout tiie paper we calculate median errors like 
l.253o-/^f(N), where <T is the standard deviation of the mean 
and N is the number of values used. We assume that the quantities 
studied here follow a normal distribution. 



not included in the following analysis. For a limit of 40 pc we 
probe the whole jet for 60% of the CJF sources (72% for BL 
Lacs, 60% for radio galaxies, and 56% for FSRQs), while for 
approximately 17% of the CJF sources more than half of the jet 
is excluded from further analysis. 

Adopting a core distance limit somewhat reduces the avail- 
able number of components and therefore possibly increases 
the scatter in the following statistical analysis. Moreover, this 
limit implicitly gives a higher weight to more strongly beamed 
sources. This bias however applies equally to both BL Lacs 
and FSRQs, and should therefore not influence the compari- 
son between the two classes. The higher number of compo- 
nents "missed" for the FSRQs because of this limit is mainly 
due to their redshift distribution and therefore using redshift bins 
should minimize the bias to our results. 

We investigate how, if at all, a jet resembles a sinusoid (as 
in the case of S5 1803+784). In order to do this we calculate 
the M.I. from the jet ridge lines of our sources. The M.I. is cal- 
culated for the whole jet (we do not apply the core separation 
limit of 40 pc), as we are interested in the morphology of the en- 
tire jet ridge line, rather than a localized property. Additionally, 
a large number of components is needed for a robust interpreta- 
tion of the M.I., we therefore take into account only epochs with 
at least three components identified. Finally, we define the max- 
imum M.I. value for each source across the available epochs. 
Fluctuations of the M.I. of a source between epochs are low. 

Under the above constraints, we find that there are in total 
46 sources with M.I. > 0.5. This translates to 37% of the total 
number of sources used here. BL Lacs show more often jets with 
sinusoid-like morphologies, compared to both FSRQs and RGs. 
More than half of the FSRQs (~ 51%) have an M.I. value of 0, 
while two thirds (~ 66%) have an M.I. value lower than 0.5. For 
RGs the respective percentages are ~ 41% and ~ 66%. In con- 
trast to this, only one third of the BL Lacs (~ 33%) have M.I.=0 
and less than half (~ 44%) of them have M.I.< 0.5. Figure 4 
shows the histogram of the M.I. for the different types of ob- 
jects. BL Lacs show their maximum in the [0.4,0.6) bin. In con- 
trast, both FSRQs and RGs have their maxima in the [0,0.2) bin. 
In Table 2 the statistical properties of the different types of ob- 
jects are given. It can be seen that the average M.I. values for BL 
Lacs and FSRQs agree within the Icr errors. When considering 
the median values for the three classes of objects, the difference 
becomes more apparent, with BL Lacs and FSRQs differing at 
> 3cr level. It should be noted that a robust quantitative treatment 
of the M.I. is problematic, given the severe restrictions imposed 
by the small number of components per jet, as well as the limited 
temporal resolution of the observations used here. 



Table 2. Statistical properties of the monotonicity index (M.I.) 
distribution for FSRQs, BL Lacs, and RGs. Average and median 
are calculated for all sources with at least three components in 
their jet. Only the maximum M.I. values are considered for each 
source. 



Monotonicity Index (max) 



Types 


FSRQ 


BL 


RG 


# 


77 


18 


29 


Average 


0.29 


0.36 


0.35 


Error 


0.03 


0.06 


0.06 


Median 





0.500 


0.330 


Error 


0.004 


0.014 


0.012 
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Fig. 4. Monotonicity Index (M.I.) distribution for BL Lacs 
(dashed line), FSRQs (solid line), and RGs (grey blocks). The 
histogram has been normaUzed to area unity. 



In this somewhat qualitative context BL Lacs appear to more 
often show sinusoid-like curved jets, as seen for S5 1803+784 
and 0716+714. 



1 . In Table 3 we show the characteristic statistical parameters for 
the CJF jet ridge line widths for all sources with available red- 
shift, as well as sources in the redshift bin [0,1]. Given the above 
analysis concerning the distance constraints, in Table 3 we give 
values also for the case where a 20 pc constraint for the maxi- 
mum distance between the 6,„ax and 0„„„ components is assumed. 

The jet ridge line width appears to be independent of red- 
shift. To test this further, we use the sub-sample of CJF FSRQs 
to calculate average and median width values for different red- 
shift bins (Fig. 5). Out to z-2.5 the jet ridge line width remains 
fairly stable, with the median value showing larger changes. For 
the last two bins, the number of sources contained drops drasti- 
cally (as reflected by the much larger error bars) and therefore it 
is diflicult to conclude whether the strong increase in both aver- 
age and median values is a true or a spurious effect. It should 
be noted that this sharp increase coincides with the epoch of 
maximum nuclear and star-formation activity in the Universe (at 
z ~ 2 - 3, e.g., Hasinger 1998; Madau et al. 1996). Focusing 
on the redshift bin [0,1], from Table 3 we see that BL Lac jet 
ridge lines appear significantly wider than their FSRQ counter- 
parts (20.4 ± 2.0°, compared to 11.1 ± 0.8°; > 4cr). The same 
behavior is seen when looking at the median values (albeit at 
the 2cr level). Considering the values when the maximum dis- 
tance limit between the 0,„ax and 6mm components is assumed, 
BL Lacs remain substantially (at a ~ 4cr difference level) wider. 
Radio galaxies show similar jet widths to BL Lacs within the Icr 
uncertainties. 



4.2. Apparent jet width, P, statistics 



In Sect. 3, we have already briefly discussed the concept of lo- 
calized width and how this relates to the measure we are going 
to investigate here. In order to better understand this relation, 
we study how the distribution of the widths of the CJF jet ridge 
lines relates to the maximum distance between the two compo- 
nents used to calculate this width. For one epoch of the source 
S5 1803+784 and one epoch of the source 2200+420, a width of 
~ 180 degrees is found. This is probably due to a problem in the 
original component fitting and we therefore exclude these two 
epochs in this analysis. 

We calculate non-zero jet widths for a total of 557 epochs 
(multiple epochs per source). Of these, 277 are larger than 10 
degrees, while 125 are larger than 20 degrees. This translates to 
22.4% of the sample possessing jets with widths exceeding 20 
degrees. 

We calculated an average jet ridge line width for our sources. 
Without imposing any constraints on the maximum distance 
between the O,,,^^ and 0,„,„ components used, we find that BL 
Lacs show significantly wider jet ridge lines (average value of 
19.3 ± 1.8°), compared to FSRQs (12.0 ± 0.4o; 4cr difference). 
RGs do not show significant differences in their mean width 
compared to BL Lacs. A Student's t-test gives a significance of 
>99.99% that BL Lacs and FSRQs show different mean values. 
After imposing a constraint of 20 pc for the maximum distance 
between the 6,nax and 0,„i„ components, the BL Lacs still show 
statistically wider jet ridge lines. Low number statistics does not 
allow us to draw a robust conclusion for more localized widths 
(at ~ 5 pc linear scales) of our CJF jet ridge lines. 

We also investigated sources in the redshift bin [0,1]. 
Although the redshift span of the CJF goes out to almost 4, of 
the 32 BL Lacs included in the CJF, only 23 have redshift infor- 
mation. Of these, 20 are below redshift 1 . Therefore for the fol- 
lowing we will focus our investigation to sources up to redshift 



-Average 
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Fig. 5. Jet width for CJF FSRQs, as a function of redshift. We 
use 0.5 redshift binning to calculate averages (continuous line) 
and median values (dashed line). Number of sources per bin: 38 
(0-0.5), 53 (0.5-1), 106 (1-1.5), 61 (1.5-2), 55 (2-2.5), 18 (2.5-3), 
12 (3-4). 



In Fig. 6 (left) we show the distribution of jet ridge line 
widths for BL Lacs, FSRQs, and RGs, in the redshift bin [0,1]. It 
can be seen that all three classes show similar distributions, with 
their maxima situated at around 10 degrees. The distribution of 
BL Lac appears to be wider and extending to larger widths, com- 
pared to FSRQs. FSRQs show a rather more contained distribu- 
tion to lower width values than BL Lacs. For the 125 epochs 
where CJF sources show apparent jet widths > 20 degrees, we 
find 14 BL Lacs, 27 FSRQs, and 13 RGs. When taking into ac- 
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Table 3. Characteristic statistical values concerning the width of the jet ridge line of sources with measured redshifts. We give 
average and median values with uncertainties, maximum, and minimum values for FSRQs, BL Lacs, and RGs. The noted distance 
pertains to the maximum distance between the two components used to calculate the jet ridge line width (see text for details). 



Jet Width (°), < z 



Types 


FSRQ 




BL 


KG 


Distance (pc) 


All 


<20 


All 


< 20 


All 


<20 


# 


343 


251 


63 


51 


95 


72 


Average 


12.0 


10.8 


19.3 


18.2 


16.3 


16.5 


Error 


0.4 


0.5 


1.8 


1.8 


1.0 


1.2 


Median 


8.7 


7.0 


12.5 


12.5 


13.0 


13.2 


Error 


0.5 


0.6 


2.2 


2.2 


0.7 


0.8 


Max 


86.4 


86.4 


93.2 


51.9 


72.8 


72.8 


Min 


0.1 


0.1 


0.7 


0.9 


0.2 


0.2 


Jet Width (°), 


< z < 1 












Types 


FSRQ 




BL 


RG 


Distance (pc) 


All 


<20 


All 


<20 


All 


<20 


# 


91 


76 


58 


49 


83 


62 


Average 


11.1 


10.5 


20.4 


18.4 


16.6 


16.4 


Error 


0.8 


0.9 


2.0 


1.9 


1.2 


1.4 


Median 


7.5 


6.5 


13.8 


12.5 


12.3 


10.8 


Error 


1.0 


1.4 


2.5 


2.4 


1.5 


1.7 


Max 


47.7 


47.7 


93.2 


51.9 


72.8 


72.8 


Min 


0.1 


0.1 


0.9 


0.9 


0.2 


0.2 



count the total number of each type of object, we calculate that 
47.3% of BL Lacs have jets wider than 20 degrees, as opposed 
to only 13.6% and 25% for FSRQs and RGs, respectively. 

We apply the two-sample Kolmogorov-Smirnoff (K-S) test 
to our data to see whether the three classes of AGN, BL Lacs, 
FSRQs, and RGs, are indeed different in their jet ridge line 
widths. Comparing the jet ridge line widths of BL Lacs and 
FSRQs (no redshift binning)the K-S test gives a 0.7 ■ 10"^% 
probability that these two sub-samples originate from the same 
parent distribution. Comparing FSRQs with RGs, we get a prob- 
ability of 0.7%. The K-S test does not reject the null hypothe- 
sis (22.1% probability) concerning the comparison between BL 
Lacs and RGs. We also compare the sub-samples of BL Lacs, 
FSRQs, and RGs in the redshift bin [0,1]. In this case the K- 
S test gives a probability of 1% that BL Lacs and FSRQs are 
drawn from the same parent sample. For BL Lacs and RGs the 
test again does not reject the null hypothesis. Finally, we also 
apply the K-S test after we assume the 20 pc distance constraint 
between the 6max/0mm components. In this case we get a 4.4% 
probability for BL Lacs and FSRQs, while for BL Lacs and RGs 
the test once more does not reject the null hypothesis. The latter 
comparison is affected by low-number statistics. We can con- 
clude that BL Lacs show a significantly different distribution of 
jet ridge line widths compared to FSRQs. In both average and 
median values, BL Lacs show substantially wider jets than both 
FSRQs and RGs. 

4.3. Apparent jet width evolution, AP, statistics 

As mentioned before, component identification is unimportant 
when working with jet ridge lines. While this is true, the ac- 
tual detection of the same number of components across differ- 
ent epochs does bear some importance for our results. The effect 
of a vanishing component (either because of diminishing flux. 



bad observing conditions, or otherwise flawed data) will induce 
a spurious width change, affecting the overall statistics of the 
sample. Although in the case of diminishing flux, one could ar- 
gue that this change is indeed intrinsic and therefore important, 
usually the latter two effects dominate. Moreover, as we already 
mentioned in Sect. 3, the given sensitivity and dynamic range 
of the observations define our results. For much more sensitive 
observations, jets would probably appear rather wider and with 
a more complex structure than what we find here. 

To take this effect into account we flag our data accordingly. 
The best quality data ("OK") are those for which a width change 
is calculated between epochs with the same number of compo- 
nents identified per source. As moderate quality data ("!") are 
flagged those for which the width change is calculated between 
epochs with + 1 number of components identified per source. If 
the number of components varies by 2 or more, the data are 
flagged as bad ("NO"). In this way we can still take into ac- 
count the ejection of new components. Within the time spans 
between most epochs (~ 2 years) the appearance of more than 
one new component is not expected. In Table 4 the number of 
epochs for each class of object with each of the corresponding 
flags is shown. As can be seen, FSRQs and RGs have the best 
quality data, with only a 24% and 27% respectively of the data 
flagged as not "OK" ("!" and "NO" flag). In contrast BL Lacs 
have almost half their data (47%) flagged accordingly. This re- 
flects the difficulty of a consistent component model-fitting of 
BL Lac jets across epochs. This is a result of (1) the variable 
nature of these objects, (2) the intrinsically fainter nature of BL 
Lacs, and (3) the possible shorter length of their jets that leads to 
confusion effects between components (compared for example 
to FSRQ jets). 

We calculate a non-zero maximum apparent width evolution 
for 180 sources (sources flagged as "NO" are not counted here). 
Of there, 5 1 show an evolution of their width larger than 4 deg/yr. 
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Fig. 6. Jet ridge line width (left) and width evolution (right) distributions for BL Lacs (dashed line), FSRQs (solid line), and RGs 
(grey blocks). For this histogram sources with redshifts in the [0,1] bin are used. Sources flagged as "NO" are excluded. 



Table 4. Three different types of flagging for our data, along with 
their definition. For each category the corresponding number for 
each class of object is shown. For an epoch i, each source has A^, 
number of components identified. 



Definition 



# of epociis 
FSRQ BL RG 



OK 
! 

NO 



Ni = Nm 
N, = Ni,, ± 1 
N, > M+i ± 2 



173 26 54 
53 21 15 

2 2 5 



while 17 have a Af„,„., > 8. We thus find that 28.3% of the 
sources show significant evolution of their apparent jet width. 

We compare the maximum values of the jet ridge line width 
evolution per source between BL Lacs, quasars, and radio galax- 
ies, using no maximum distance limit between the 0,„„ i and 0,„i„ 
components for the width calculation and excluding data flagged 
as "NO". The differences in average and median values for 
the different classes are not found to be statistically significant 
(< 3cr). BL Lacs do however show a trend of having stronger 
width evolution than the other two classes. When checking the 
maximum and minimum values of each sub-sample, it becomes 
obvious that FSRQs and BL Lacs show much more extended 
distributions compared to RGs. 

In Fig. 6 (right) we show the distributions for the maxi- 
mum jet ridge line width evolution for FSRQs, BL Lacs, and 
RGs. We use a redshift bin of [0,1] and exclude all data flagged 
as "NO" (this corresponds to the last three columns of the left 
part of Table 5). BL Lac objects show a wide distribution of jet 
width evolution values, with a primary maximum at 3 degrees/yr 
and extending out to 18 deg/yr Radio galaxies show a some- 
what similar behavior with a plateau between 2 and 4 degrees/yr 
Quasars, in contrast, show a confined distribution with a strong 
maximum around 3 degrees/yr and all values contained between 
and 10 degrees/yr. We apply the K-S test to see whether 
these distribution are significantly different. Concentrating on 
the comparison in the [0,1] redshift bin, the K-S test does not 
reject the null hypothesis (33.8% probability). Calculating the 



relative occurrence of the different classes, we find that 12.5% 
of BL Lacs show t^Pmax > 8, as opposed to 4% and 5.7% for 
FSRQs and RGs, respectively. 

4.4. Apparent jet linear evolution, At, statistics 

Focusing on the final measure of the jet ridge Une kinematics, as 
studied in this paper, we want to investigate whether the station- 
arity of components, as observed in the case of S5 1803-1-784, is 
commonplace among other BL Lacs. We use the linear evolution 
measure, as described in Sect. 3, to do this. We follow the same 
procedure as previously, to check the statistics of the individual 
classes. 

In Table 5 (right) we give the statistical properties of the total 
linear evolution of the jet ridge line (measured in parsecs per 
unit time and component) distributions for FSRQs, BL Lacs, and 
RGs. On average, both BL Lacs and RGs show weaker evolution 
of their jet ridge lines compared to FSRQs (~ 2cr difference). 
Looking at the median values, BL Lacs and RGs show the least 
evolution compared to FSRQs (0.24 + 0.06, 0.28 + 0.04, and 
0.414+0.020 pc/yr/comp respectively; ~ 3cr difference). We also 
study the statistics of the sub-sample of CJF sources in the [0,1] 
redshift bin (also in Table 5). The behavior remains the same as 
before, with differences between BL Lacs and FSRQs becoming 
more pronounced (with a ~ 4cr difference in median values). 

We now turn to the actual distribution of the total jet linear 
evolution for the three types of objects, shown in Fig. 7. We once 
again only take into account sources in the redshift bin [0,1]. 
Quasars show a pronounced maximum around 0.63 pc/yr/comp, 
extending out to 2 pc/yr/comp. In contrast, both BL Lacs and 
RGs show their distribution maxima around 0.25 pc/yr/comp. 
Quasars appear to have a wider distribution of total linear evo- 
lution values, also showing the highest maximum values among 
all three types (see Table 5). Interestingly, RGs and BL Lacs, 
in total, show very similar distributions. Once again we employ 
the K-S test to compare the distributions. For the whole sample 
(independent of redshift constraints), the K-S test gives a proba- 
bility of 4.3% that BL Lacs and FSRQs are drawn from the same 
parent population. For BL Lacs and RGs the test does not reject 
the null hypothesis. Focusing on the < z < 1 sub-sample, we 
get a marginal 5.4% probability that BL Lacs and FSRQs stem 
from the same parent sample. These results give a positive an- 
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Table 5. Statistical properties of the maximum jet width evolution (left) and jet hnear evolution (right) distributions for FSRQs, BL 
Lacs, and RGs. Average, median , maximum, and minimum values are calculated for both all sources, as well as for sources in the 
redshift bin [0,1]. Sources flagged as "NO" are excluded. Only the inner part of the jet is considered (< 40 pc). 



Jet width evolution ( 


max) (deg/yr) 




Jet linear evolution (pc/yr/comp) 








All 






0<z< 1 






< z 






< z < 1 




Types 


FSRQ 


BL 


RG 


FSRQ 


BL 


RG 


FSRQ 


BL 


RG 


FSRQ 


BL 


RG 


# 


114 


26 


35 


29 


19 


26 


171 


25 


41 


44 


21 


35 


Average 


3.44 


3.9 


3.3 


3.5 


4.1 


3.7 


0.470 


0.38 


0.34 


0.51 


0.38 


0.32 


Error 


0.23 


0.6 


0.3 


0.4 


0.7 


04 


0.016 


0.05 


0.03 


0.04 


0.05 


0.03 


Median 


2.26 


2.3 


2.5 


2.42 


2.8 


2.73 


0.414 


0.24 


0.28 


0.50 


0.24 


0.25 


Error 


0.29 


0.7 


0.4 


0.5 


0.9 


0.5 


0.020 


0.06 


0.04 


0.04 


0.06 


0.04 


Max 


27.11 


18.32 


11.84 


24.39 


18.32 


11.84 


1.910 


1.081 


0.989 


1.910 


1.081 


0.989 


Min 


0.16 


0.34 


0.14 


0.41 


0.12 


0.51 


0.028 


0.038 


0.028 


0.028 


0.038 


0.028 



swer as to whether BL Lacs show less absolute linear evolution 
of their jets compared to their FSRQ counterparts. A somewhat 
surprising result concerns the comparison between BL Lacs and 
RGs. 

4.5. Jet ridge line properties and variability 

Karouzos et al. (2010) gathered variability information from the 
literature for the CJF sources, focusing on the long timescale 
(~ 1 - 10' years) variability, often linked to a helical jet structure 
and/or a binary AGN core. We use this information to define a 
subsample of 40 CJF variable sources and study their properties 
in terms of the measures discussed in Sect. 3. Given the small 
number of sources/epochs for this investigation, we impose no 
additional constraints to our sub-sample (e.g., redshift bins). 

The sub-sample of variable sources exhibits both higher av- 
erage (18.3± 1.3°, compared to 14.4 + 0.4°; ~ 3cr difference) and 
median (15.7 ± 0.9°, compared to 9.70 ± 0.26°; 6cr difference) 
width values. In Fig. 8 (left) we plot the jet ridge line width dis- 
tribution for the sub-sample of variable sources, as well as for 
the rest of the CJF. We see that the parent (non-variable^) CJF 
sample shows strong maximum in the lowest bin [0°, 8°] with a 
decreasing trend. In contrast, the sub-sample of variable sources 
show two maxima in bins [0°,8°] (main) and [16°, 24°] (sec- 
ondary), with a secondary to primary ratio of ~ 0.9. Variable 
CJF sources indeed show a more extended distribution of jet 
ridge line widths, shifted towards higher values compared to the 
parent CJF sample. A K-S test gives a very low (0.3%) proba- 
bility that the two sub-samples are drawn from the same parent 
population. 

BL Lacs appear to maintain the same average and median 
width as their parent sample (within errors). FSRQs however ap- 
pear to have both much higher average (22.0 ± 2.7°, compared 
to 10.8 + 0.5°; 4cr difference) and median values than their par- 
ent sample (17.8 ± 1.6°, compared to 7.00 ± 0.28°; 4cr differ- 
ence). Variable RGs showmarginal evidence for lower average 
(12.9 + 2.1°, compared to 16.5 ± 1.2°) but similar median value. 

For AP, compared to the parent sample, the sub-sample of 
variable sources exhibits both higher average (5.5+0.9° /yr, com- 
pared to 3.51 + 0.18° /yr) and median values (3.1 +0.4° /yr, com- 

' We underline the fact that the number of variable sources is explic- 
itly regulated by the availability, or lack thereof, of the appropriate ob- 
servations (long-timescale variability can be uncovered only by means 
of extensive monitoring of a source). Some of the sources labeled here 
as "non-variable" are probably variable but not observed as such. 



pared to 2.37 + 0.11° /yr), albeit at a low (2cr) significance level. 
Variable BL Lacs appear to maintain both the same average 
width as their parent sample (within errors), as well as the same 
median value (2.5 ± 0.5° /yr, compared to 2.3 + 0.3° /yr). FSRQs 
appear to have both higher average (8.4 ± 2.3° /yr, compared to 
3. 44 + 0.22° /yr; 2cr diff'erence) and median values (5.0± 1.0°/^^ 
compared to 2.26 + 0.13°/3'r; 2.5sigma) than their parent sam- 
ple. Variable RGs show lower average (1.9 ± 0.6°/3'r, compared 
to 3.3 ± 0.3°/3'r) but similar median values (1.6 ± 1.0°/3'r, com- 
pared to 2.48 ±0.20° /yr). 

For M, compared to the parent sample, the sub-sample of 
variable sources does not exhibit significantly different average 
and median values. For both samples the total maximum of the 
distribution occurs in the [0.2,0.4] bin, with the variable CJF 
sources however showing a less strong maximum compared to 
the non-variable sample (see Fig. 8, right). In contrast to the non- 
variable sources, the variable sample shows a secondary maxi- 
mum in the bin [0.8,1] pc/yr/comp implying a broader distribu- 
tion of M for variable sources. The K-S test does not reject the 
null hypothesis (23.8% probability). BL Lac objects appear to 
maintain the same average and median linear evolution as their 
parent sample (within errors). Variable quasars appear to have 
both higher average (0.615 + 0.063, compared to 0.470 ± 0.016; 
> Act difference) and median values (0.53 ± 0.04 pc/yr/comp, 
compared to 0.414 + 0.020 pc/yr/comp) than their parent sam- 
ple. Variable RGs show similar average values (within errors) 
but lower median value (0.135 + 0.003 pc/yr/comp, compared to 
0.28 ± 0.04 pc/yr/comp). 

There seems to be a strong link between the kinematic prop- 
erties studied here and the variability observed in these sources. 
For all three measures BL Lacs appear to maintain their statisti- 
cal behavior between variable and parent sample, whereas vari- 
able quasars tend to lean towards a "BL Lac-like" behavior It 
should be noted however that, especially for RGs, this compari- 
son is problematic due to the small number of sources (3 RGs). 

4.6. Jet kinematics and source luminosity 

In addition to the above, we also look for a possible correlation 
between the kinematic and morphological measures described in 
Sec. 3 and the radio luminosity of the sources, a value probably 
less heavily dependent on the jet structure/morphology. 

In Fig. 9 we show the total linear evolution of the CJF jet 
ridge lines as a function of the core luminosity of the source at 5 
GHz (from VLA observations). It can be immediately seen that a 
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Fig. 7. Distribution of the total linear evolution of the CJF jet ridge lines for FSRQs (solid line), BL Lacs (dashed line), and RGs 
(grey blocks). We use sources with redshift lower than 1. 



Table 6. Statistical properties of the jet width, maximum jet width evolution, and jet linear evolution for the sub-sample of variable 
CJF sources (see text for definition of variability). Average and median are calculated for all sources, independent of redshift. 
Sources flagged as "NO" are excluded. Only the inner part of the jet is considered (< 40 pc). For each measure the upper row gives 
average values, while the median values are found in the lower rows. 





All 


FSRQ 


BL Lac 


RG 


p 


18.3±1.3 


22.0±2.7 


17.1±1.7 


12.9±2.1 


(deg) 


15.7±1.6 


17.8±3.3 


13.5±2.1 


13.4±2.6 


AP 


5.5±0.9 


8.4±2.3 


4.5±0.9 


1.9±0.6 


(deg/yr) 


3.1±1.1 


5.0±2.9 


2.5+1.1 


1.6±0.7 


M 


0.49±0.05 


0.61±0.06 


0.36±0.07 


0.35±0.17 


(pc/yr/comp) 


0.38±0.06 


0.53±0.08 


0.23±0.08 


0.14±0.21 



correlation exists between radio luminosity and total linear evo- 
lution of the jet ridge line. This is seen as an upper envelope, 
above which no sources are found. More strongly evolving jets 
inhabit brighter sources. Conversely, low luminosity sources ap- 
pear to show the smallest values of linear displacement of their 
jets. A similar correlation has been reported between maximum 
apparent jet-component speed and radio luminosity (Cohen et al. 
2007; Britzen et al. 2008; Lister et al. 2009). Britzen et al. (2008) 
discuss this correlation between apparent speeds and luminosity 
in terms of a redshift induced effect. As can be seen in Fig. 9, 
there is an indication that this is not an effect caused by redshift 
(e.g., CJF sources at the highest redshifts (> 3), albeit only being 
a few, all appear to have rather weak linear evolution of their jet 
ridge lines). Further arguments against a redshift induced effect 
dominating this correlation can be found in Britzen et al. (2008). 
It should be noted that small-number statistics might be affect- 
ing our results. Both for low luminosity sources, as well as high 
redshift sources, we have a small number of objects available. 
That hinders the global validity of what is seen in Fig. 9. 

Concerning differences between FSRQs, BL Lacs, and RGs, 
no strong effects can be seen in this context. The low-luminosity 
tail of the distribution is predominantly populated by RGs and 
BL Lacs, with FSRQs mainly dominating the higher luminosi- 
ties. Cohen et al. (2007) use Monte Carlo simulations to show 
that such an envelope can be explained in terms of "aspect 
curves", essentially supporting the beaming model of relativistic 
jets as an explanation for this effect. In Fig. 9 (right) we plot the 



same correlation but only for the RGs in our sample. They form 
the same upper envelope, with more strongly evolving jets ap- 
pearing in more luminous RGs. In the context of the unification 
scheme, where RGs are believed to be observed at large angles 
to their jet axes. Fig. 9 indicates that the correlation between to- 
tal linear evolution of the jet ridge line and the luminosity can 
not be attributed only to geometric effects. 

5. Discussion 

Given the flux-limited nature of the CJF sample, we are inter- 
ested in investigating possible factors influencing the results pre- 
sented in this paper. These can be briefly summarized in the fol- 
lowing: 

- relative number of FSRQs and BL Lacs 

- CJF redshift range 

- jet component identification 

- projection effects 

We shall address each of these points separately. 

5.0.1. Total and absolute source numbers 

Compared to other samples of similar nature (e.g., MOJAVE), 
the CJF contains a larger number of sources, enabling the study 
of individual populations or sub-groups of sources (i.e., FSRQs, 
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Fig. 8. Jet ridge line width (left) and linear evolution (right) distributions for the sub-sample of variable sources from Karouzos et al. 
(2010) (black line) and for the rest (non-variable) CJF sample (grey blocks). Due to the big difference in absolute numbers for the 
two samples, we have normaUzed to unity surface area. No redshift or distance constraints are assumed. 




Fig. 9. Total linear evolution of the jet ridge lines for the BL Lac and RG sub-samples (up; left and right, respectively) and the FSRQ 
sub-sample (down) as a function of the total source luminosity at 5 GHz. Quasars are shown with squares, BL Lacs with circles, 
and RGs with triangles. We use a color gradient to denote different redshifts. 



BL Lacs, etc.). It should be noted that, given the selection crite- 
ria of the sample and of course the relative cosmic abundance of 
these sources, quasars dominate the sample, with only a few BL 
Lacs and radio galaxies in comparison. This introduces a cer- 
tain degree of small number statistics uncertainties for some of 
the comparisons in this paper. Future radio surveys with next- 
generation instruments will surely allow the definition of much 
larger samples that will address this problem. 



5.0.2. Redshift distribution 

The large redshift span of the CJF, as well as the uneven sam- 
pling of that range with respect to the three types of AGN in- 
fluence our results. BL Lac objects are over-represented in the 
redshift range to 2. Radio galaxies follow a similar behav- 
ior, with most sources in the redshift range [0,1], but a few 
sources at higher redshifts. Contrary to that, FSRQs show are 
over-represented in the redshift bin [1,2] and extending to the 
highest redshift values of the sample. 
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As we briefly discussed before, the fixed sensitivity and res- 
olution of the arrays used allows us to have a direct comparison 
between the different CJF sources. On the same time however, it 
also essentially leads us to observe different scales of the jet for 
objects at different redshifts. It is expected that the properties, 
kinematics, and morphology of a jet are strong functions of the 
separation from the core. We therefore adopt redshift bins for all 
of the comparisons we do. Given the distribution of redshifts for 
the CJF the only bin we can use without losing most of the BL 
Lacs is the [0,1] one"*. 

5.0.3. Component identification 

The way that we study the width and width evolution of the jet is 
completely independent of component identification. The linear 
evolution of the jet ridge line on the other hand is an excep- 
tion to this. For this measure we are forced to follow the given 
identification of components (from Britzen et al. 2007b) as we 
need some reference point to calculate absolute displacements. 
Despite the inevitable coupling of prior component identifica- 
tion to our A/" measure, the fact that we include in this measure 
all components and all epochs ensures a treatment of the jet as 
a whole, evening out peculiarities or possible misidentifications 
of individual components. 

Our measures are however sensitive to the number of compo- 
nents identified at a certain epoch. This is especially relevant for 
BL Lacs, given the nature of their jets. As we discussed shortly 
previously, BL Lacs are extremely variable sources. In parts this 
owns to the presumably small angle to our line of sight, resulting 
to beaming effects that change the flux of both the core and the 
jet. Secondly, we have shown that BL Lacs show indications of 
shorter jets than both FSRQs and RGs. This translates to higher 
difficulty in consistently model-fitting the jets of BL Lacs, as 
blending effects become more prominent as the viewing angle 
decreases. Moreover, given the strong variability of both core 
and jet, components might simply vanish as a result of dimin- 
ishing luminosity or insufficient dynamic range of the observa- 
tions. A further way to demonstrate this is by using the quality 
classification for jet components used by Britzen et al. (2008). 
Jet component proper motions are classified as Ql for the best 
quality data, and Q2 and Q3 for diminishing quality. The ratio 
N{Ql)/{N(Q2) + N{Q3)), where N(x) is the number of x quality 
components of that type of source, is 0.63 for FSRQs, 0.80 for 
RGs, and 0.41 for BL Lacs. Returning to our original point, we 
see therefore that the number of components across epochs in 
BL Lacs jets is more variable than both FSRQs and RGs. This 
introduces some uncertainty to our results. 

5.0.4. Projection effects 

The fact that the values calculated here are all projected onto the 
plane of the sky introduces some uncertainty and consequently 
a scatter in this statistical investigation. Although for this analy- 
sis and for the following discussion we have adopted the general 
AGN unification scheme, where BL Lacs and FSRQs are at the 
smallest viewing angles and radio galaxies are at larger view- 
ing angles, the dispersion of the actual angles within each class 
introduces the aforementioned scatter in our statistics. 



* It sliould also be noted that given the featureless spectrum, for 
which BL Lacs are famous, leads to more than 34% of the CJF BL 
Lacs to not have available redshift information. When using redshift 
bins these objects are obviously excluded. 



An alternative path would be to try and deproject the jet val- 
ues calculated here. Although certainly possible for a number 
of the CJF sources where a Doppler factor has been calculated 
(e.g., Britzen et al. 2007a, it would introduce further degrees of 
uncertainty to our analysis that are not so easily constrained. As 
we are mainly interested in the comparison between FSRQs and 
BL Lacs rather than the calculation of intrinsic jet properties, the 
deprojection is not essential for the results of this paper 

5.1. Possible explanations 

It is of great interest to try and explain the kinematic behavior 
seen in a large number of the CJF (combining the percentages 
calculated from both apparent width and apparent width evolu- 
tion we get an approximately 30% of the sample showing wide 
jets that change their width strongly) and the apparent preference 
for BL Lacs to show this behavior over FSRQs. One obvious fac- 
tor that should play a deciding role for the kinematics of jet com- 
ponents is the viewing angle under which a source is observed. It 
is known that BL Lacs and FSRQs are seen jet-on, at the smallest 
viewing angles, with steeper-spectrum quasars and radio galax- 
ies being sources observed at progressively larger viewing an- 
gles (e.g., Antonucci 1993; Urry & Padovani 1995). Could a 
viewing angle difference explain the kinematic differences seen 
in some of the CJF sources and in particular observed between 
BL Lacs and FSRQs? Assuming that jet components follow bal- 
listic, linear paths, one can expect that viewed at smaller view- 
ing angles (smaller than the critical angle 1 /y, with y being the 
Lorentz factor of the flow), the components are observed to cover 
smaller distances, than if seen edge-on. Therefore, the slower 
speeds, and hence smaller total linear evolution of their jet ridge 
lines, observed for BL Lacs could be explained in terms of a 
systematically smaller viewing angle for BL Lacs, compared to 
FSRQs. 

However, Hovatta et al. (2009) (based on an investigation 
initiated by Lahteenmaki & Valtaoja 1999) used long timescale 
variability of AGN, together with jet kinematics, to derive 
Doppler factors and consequently viewing angles of a sample 
of 87 AGN. They found that FSRQs actually show smaller 
mean viewing angles than BL Lacs (a result also found by 
Lahteenmaki & Valtaoja 1999). Although there are certain limi- 
tations to the method used by the authors, their results are diffi- 
cult to reconcile with a viewing-angle-dependent explanation of 
our results. 

Our results would therefore indicate the need for an addi- 
tional, potentially geometric, effect in play. This would be in 
agreement with Cohen et al. (2007), who postulate that low- 
speed components in FSRQ and BL Lac jets appear so because 
their pattern Lorentz factor is lower than the bulk Lorentz factor 
of the jet. Alternatively, it should be considered whether there is 
some systematic bias in the way Doppler factors are estimated, 
that would lead to an over- or underestimation of the viewing 
angles of BL Lacs and FSRQs, respectively. 

The large jet widths and jet widths changes are more diffi- 
cult to explain. An additional mechanism or effect must be in- 
troduced to produce a wide jet. Precession of the jet axis, or as- 
suming that the components follow non-ballistic and non-linear 
paths can lead to such jet properties (e.g., Steffen et al. 1995; 
Gong 2008; Roland et al. 2008; Gong et al. 2011). In these 
models, jet components follow highly curved trajectories, which, 
viewed face-on, would give the impression of a wider distribu- 
tion of jet component position angles. A precessing jet, or rather 
a precessing jet nozzle, would imply that different components 
follow different trajectories. That would result in a changing 
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jet component position angle distribution, more pronounced at 
smaller viewing angles. That more BL Lacs show M.I. closer to 
one compared to FSRQs, implying a highly curved, sinusoid-like 
ridge-line, lends additional support to this scenario. Following 
this train of thought would then lead us to the conclusion that ei- 
ther (1) there is a sub-set of FSRQs and BL Lacs (~ 30% in our 
sample) for which non-ballistic and/or precession effects play an 
important role, or more generally that (2) for these sources the 
viewing angles are below some critical angle 6c that allow us to 
witness and thus characterize the helical or "non-ballistic" struc- 
ture common in all radio-AGN jets, but which would otherwise 
be blended out by the projection effects at viewing angles > Oc- 

In a following paper, we shall use an expanded version of the 
helical jet model of Steffen et al. (1995) to evaluate the statistical 
results presented here, in terms of the viewing angle (and hence 
beaming) effect, as well as a possible helical geometry of the jet. 

A final scenario that needs to be discussed is whether a sys- 
tematic difference in the Lorentz factors y between FSRQs and 
BL Lacs could produce the effects observed here. It has been 
argued that BL Lacs show smaller Lorentz factors than FSRQs 
(e.g., Morganti et al. 1995; Urry & Padovani 1995; Hovatta et al. 
2009). If this is true (and not a selection bias in the samples used) 
then that would naturally explain the slower components in BL 
Lacs. The above in turn ties in to the currently accepted unifi- 
cation scheme (e.g., Urry & Padovani 1995). In that paradigm, 
BL Lacs and FSRQs are drawn from two, presumably different, 
parent samples of Fanaroff-Riley I (low luminosity) and II (high 
luminosity) galaxies (FR; Fanaroff & Riley 1974; Padovani & 
Urry 1991; Capetti & Celotti 1999; Xu et al. 2009), respectively. 
If that is true, then the differences seen in the jet ridge line prop- 
erties of BL Lacs and FSRQs should then translate to differences 
between FRJ and FRII jet kinematics. Interestingly, studies of 
FRI and FRII kinematics have shown that both types have simi- 
lar parcec-scale jet Lorentz factors (e.g., Giovannini et al. 2001), 
contradictory to what we find here, as well as in other studies 
mentioned above, concerning the parsec-scale jet speeds in BL 
Lacs and FSRQs. 

It should be noted however that FRI sources show disrupted 
jets much earlier (i.e., closer to the core) than FRIIs (e.g., see 
Laing 1996 and references therein). Perucho et al. (2010) explain 
this in terms of decollimation through the growth of non-linear 
Kelvin-Helmholtz instability modes in FRI jets. In this context, 
the wider jets seen in BL Lacs could indicate larger amplitude 
instabilities (closer to the non-linear regime) in the jet that, given 
the lower-speed flows observed in these objects, would lead to an 
earlier disruption of the jet at kpc scales (as expected for FRIs). 

It is difficult to offer a single answer concerning an interpre- 
tation of our results. It is very likely that both viewing angle and 
Lorentz factor effects, as well as flow instabilities that depend on 
the latter, influence the appearance of BL Lac and FSRQ jets. A 
comparison between the jet ridge lines of FRIs and FRIIs to BL 
Lacs and FSRQs, respectively, would offer further insight on the 
problem. A much larger sample than the CJF would be required 
for such a study. 

5.2. Comparison with previous studies 

It should be noted that this is the first time the jet ridge lines of 
a sample of AGN this size have been explicitly and exclusively 
studied in a statistical manner. While the ridge line of a jet is not 
a new concept, so far only the jet ridge lines of individual sources 
(both galactic and stellar) have been studied (e.g., Condon & 
Mitchell 1984; Steffen et al. 1995; Lister et al. 2003; Lobanov 
et al. 2006; Britzen et al. 2009; Perucho et al. 2009; Lister et al. 



2009; Britzen et al. 2010a). The definition of the jet ridge Une 
can differ from study to study. We define the ridge line of a jet 
at a certain epoch as the line that linearly connects the projected 
positions of all components at that epoch. Alternatively, one can 
use an algorithm to more directly extract the jet ridge line from 
the VLBI map of each epoch (e.g., Pushkarev et al. 2009). One 
should however be careful of over-sampling any given map, in 
effect extracting information that is not actually there. In this 
sense, the method followed here, although perhaps simpler, is 
more robust in the quality of information used. 

The one obvious direct comparison that can be drawn is to 
the MOJAVE / 2cm sample (Lister & Homan 2005; Paper I of 
a series of 6 papers). As the work presented here is prototype 
in its conception, it is difficult to draw direct comparisons to 
the MOJAVE sample. It is of interest to compare the apparent 
velocity distributions for the MOJAVE and the CJF, as the ap- 
parent velocities are in a sense reflected in our measure of the 
total jet ridge line linear evolution. The CJF sources appear to 
be significantly slower than what is seen in the MOJAVE sample 
(e.g., compare Fig. 7 from Britzen et al. 2008 and Fig. 7 from 
Lister et al. 2009). In the distribution of apparent speeds, for the 
CJF there is a maximum at around 4c and then a turn down, 
with very few sources above 15c. In contrast, the distribution 
for the MOJAVE sources shows a maximum at 10c with a fair 
number of sources showing speeds greater than 20c. This can 
be explained in terms of the higher temporal resolution of the 
MOJAVE sample (owning to the larger number of epochs per 
source) that might allow the detection of such faster motions. 
It is interesting to note that no clear distinction has been made 
in the MOJAVE sample between BL Lacs and FSRQs, proba- 
bly due to the relatively small number of BL Lacs. Lister et al. 
(2009) however do mention that BL Lacs are more often found 
to exhibit "low-pattern speeds", i.e., components with slow mo- 
tions, significantly smaller than others in the same jet. Similar 
evidence has been previously found in different VLBI samples 
(e.g., Jorstad et al. 2001; Kellermann et al. 2004). Combined 
with the argument that BL Lac jets exhibit lower Lorentz fac- 
tors, as a result of a possible correlation between intrinsic AGN 
luminosity and jet Lorentz factor (e.g., Morganti et al. 1995), 
the above are in agreement with CJF BL Lacs showing on aver- 
age weaker linear evolution than FSRQs. It should be noted that 
the agreement with previous studies concerning the outward mo- 
tions of BL Lacs objects supports the robustness and reliability 
of our method. 

As was previously mentioned, a correlation between appar- 
ent speeds and luminosity has been studied for the MOJAVE 
sample too (Cohen et al. 2007; Lister et al. 2009). A similar cor- 
relation was observed also for the CJF sample, as was shown in 
Britzen et al. (2008). It should be noted that for both samples, 
these correlations were between luminosity and either the max- 
imum apparent speed per source (in the case of MOJAVE), or 
only for the best quality (Ql) data (for the CJF). In this paper we 
presented a similar correlation between 5 GHz core luminosity 
(as derived from VLA measurements) and the total linear evolu- 
tion of the jet ridge line. In this sense, we include speed infor- 
mation for all components identified across-epochs and therefore 
treat the jet as a whole, rather than singling out the fastest com- 
ponent. The result is however remarkably similar to what was 
previously found. Cohen et al. (2007) and Lister et al. (2009) 
discuss possible interpretations of the upper-envelope that is ob- 
served in this correlation, in terms of the beaming model. They 
conclude that there is no apparent argument as to why such an 
envelope should exist. They instead offer an alternative interpre- 
tation in terms of a link between the energy output of the AGN 
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and the intrinsic Lorentz factor of a jet. They conclude however 
that the shallow flux limit of the MOJAVE sample does not allow 
for a robust conclusion. The deeper flux-limit of the CJF (and 
the resulting more than double number of sources) allow us to 
re-address the same issue. The fact that this correlation remains 
in the CJF sample, both in the sense of a maximum apparent jet 
speed, as well as that of the linear evolution measure (in essence 
a mean jet speed), argues in favor of a link between the intrinsic 
luminosity and Lorentz factor. A definite answer would come by 
the inclusion of low -luminosity sources in such an investigation. 
It would then become obvious whether this envelope holds and 
potentially differentiate between a beaming model signature and 
that of an intrinsic connection of the energy output and the jet 
kinematics. 

A final note concerning the comparison between our results 
and these of the MOJAVE sample pertains to the observations 
frequency of each sample. Unlike the CJF, the MOJAVE uses 
VLBI observations at 15 GHz, meaning that a potentially differ- 
ent regime of the jet is probed, compared to what is studied in 
the CJF. Two-zone models, such as the "two-fluid" model (e.g., 
Sol et al. 1989, Pelletier & Roland 1989), predict a fast jet spine, 
embedded in a slower moving sheath. This would imply that the 
MOJAVE potentially probes deeper, thus faster, jet layers com- 
pared to what is seen in the CJF sample. Moreover, higher ob- 
servation frequency also translates to higher angular resolution 
and therefore to probing effects closer to the core and also en- 
abling the detection of faster moving components. In this case, 
it would be interesting to apply the methods used in this paper 
for the MOJAVE dataset, as it is expected that the curvature of 
the jet (and therefore potentially the angular properties and evo- 
lution of the jet ridge line) should play an increasingly important 
role at closer core separations. 



6. Conclusions 

We summarize our results: 

- we develop a number of tools to investigate both the mor- 
phology and kinematics of the CJF jet ridge lines: 

1. Monotonicity Index 

2. Apparent jet width 

3. Apparent jet width evolution 

4. Apparent jet linear evolution 

- Using the M.I., we find that BL Lacs jet ridge lines more 
often resemble a sinusoidal curve compared to both FSRQs 
and RGs. In contrast, 2/3 of the FSRQs have an M.l. lower 
than 0.5 (indicating fairly mono tonic jets). 

- 22.4% of the CJF sources have apparent jet widths larger 
than 20 degrees. 47.3% of BL Lacs, over 13.6% and 25% for 
FSRQs and RGs, respectively, have dP > 20 degrees. 

- BL Lacs exhibit substantially apparently wider jets. This 
effect persists under several constraints (e.g., redshift bins, 
core separation limits). Quasars show the least wide jets 
among the three classes. This supports the efi'ects seen in 
individual BL Lac objects (i.e., 1803+784, 0716+714, etc.). 

- No indication is found for a change of apparent jet width 
with redshift, at least out to z ~ 2. 

- The distribution of apparent jet ridge line widths for BL Lacs 
appears to extend towards higher values, with FSRQs mainly 
contained at lower values. Radio galaxies show a somewhat 
wider distribution of widths than FSRQs. A K-S test indi- 
cates a 0.7 ■ 10"^% probability that FSRQs and BL Lacs are 
drawn from a single parent population. 



- 28.3% of the CJF sources exhibit apparent jet width evolu- 
tion larger than 4 deg/yr. 12.5% of BL Lacs have hP,„ax > 4, 
as compared to 4% and 5.7% for FSRQs and RGs, respec- 
tively. 

- When taking into account the maximum values of appar- 
ent width evolution for the CJF sources, BL Lacs are found 
to show marginally stronger evolution of their apparent jet 
widths, with FSRQs and RGs showing similar values (within 
errors). This behavior holds for a smaller redshift bin. 

- BL Lac objects, on average, show weaker apparent lin- 
ear evolution of their jet ridge lines compared to FSRQs. 
Interestingly, RGs show even weaker linear evolution. BL 
Lac objects and radio galaxies show apparent slower moving 
jets than quasars. 

- Variable CJF source are, on average, found to exhibit appar- 
ent wider jets than the non-variable ones. A K-S test gives 
a 1 % probability that these two sub-samples stem from the 
same parent distribution. 

- Variable sources in the CJF show, on average, stronger ap- 
parent angular evolution than the non-variable ones. This is 
reflected in the fairly wider distribution of apparent jet width 
evolution angles compared to the parent (non-variable) sam- 
ple. 

- Variable sources exhibit stronger apparent linear evolution. 
This is seen also in the distributions of the two sub-samples. 

- Variable FSRQs show more "BL Lac -like" behavior with 
wider jets that also change their width more strongly. 

- A correlation, in the form of an upper envelope, is found 
between total apparent linear evolution of the CJF jet ridge 
lines and their radio (VLA, 5 GHz) core luminosities. 
Although partly attributed to relativistic effects, there are in- 
dications of a further, intrinsic effect additional to the geo- 
metric one. 

In conclusion, by statistically analyzing the CJF sample, 
we provide detailed insight concerning the morphology and 
evolution of AGN jet ridge lines. The statistical investigation 
of the CJF sources lends independent support to the different 
kinematic scenario recently seen in a number of BL Lac ob- 
jects (1803+784, Britzen et al. 2010a; 0735+178, Britzen et al. 
2010b; etc.). 25%-30% of the CJF sample show apparently con- 
siderably wide jets and strong apparent width evolution. BL 
Lac objects appear to deviate the strongest from the kinematic 
paradigm, widely accepted for blazars, of outward superlumi- 
nally moving jet components. BL Lacs appear to evolve their jet 
ridge lines (with respect to the core) less than the other source 
classes, hence indicating a slower apparent flow in their jets. On 
the other hand, they show significantly apparent wider jets, with 
a trend to change their widths more strongly, than both FSRQs 
and RGs. All these three effects have been already observed 
in the jets of 1803+784, 0716+714, and 0735+178. Viewing 
BL Lacs at small angles to their jet axis possibly allows us to 
uncover this peculiar kinematic behavior that would otherwise 
be inaccessible at larger viewing angles. However, the picture 
does not remain "simple". By studying a sub-sample of CJF 
sources identified as variable, we find that these sources actually 
show wider, more strongly evolving jet ridge lines, with variable 
FSRQs exhibiting a more "BL Lac like" behavior All the above, 
combined with the significant number of CJF sources showing 
evidence supporting a different kinematic scheme (independent 
of their classification), imply a rather universal effect, rather than 
something unique to BL Lacs. 

It should therefore be noted that the above results underline 
the fact that the notion of linear, ballistic trajectories for AGN 
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jet components usually employed until recently is a very crude 
approximation and, more often than not, deviates grossly from 
the reality. It is of great interest to uncover the process that leads 
to the properties of AGN jet ridge lines as studied in this paper. 
In a following paper, we shall use a simple jet model to address 
this point. 

Acknowledgements. M. Karouzos was supported for this research through a 
stipend from the International Max Planck Research School (IMPRS) for 
Astronomy and Astrophysics. M.K. wants to thank Lars Fuhrmann and Mar 
Mezcua, as well as an anonymous referee for insightful discussions and com- 
ments that greatly improved this manuscript. This research has made use of 
the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. This research has made use 
of NASA's Astrophysics Data System Bibliographic Services. The Very Long 
Baseline Array is operated by the USA National Radio Astronomy Observatory, 
which is a facility of the USA National Science Foundation operated under co- 
operative agreement by Associated Universities, Inc. 



References 

Agudo, I., Gomez, J., Marti, J., et al. 2001, ApJ, 549, L183 

Agudo, I., Gomez, J. L., Gabuzda, D. C, et al. 2006, A&A, 453, 477 

Antonucci, R. 1993, ARA&A, 31, 473 

Britzen, S., Brinkmann, W., Campbell, R. M., et al. 2()07a, A&A, 476, 759 
Britzen, S., Kam, V. A., Witzel, A., et al. 2009, A&A, 508, 1205 
Britzen, S., Kudryavtseva, N. A., Witzel, A., et al. 2010a, A&A, 511, A57+ 
Britzen, S., Vermeulen, R. C, Campbell, R. M., et al. 2008, A&A, 484, 1 19 
Britzen, S., Vermeulen, R. C, Taylor, G. B., et al. 2007b, A&A, 472, 763 
Britzen, S., Vermeulen, R. C, Taylor, G. B., et al. 1999, in Astronomical Society 

of the Pacific Conference Series, Vol. 159, BL Lac Phenomenon, ed. L. O. 

Takalo & A. Sillanpaa, 431-+ 
Britzen, S., Witzel, A., Gong, B. R, et al. 2010b, ArXiy e-prmts 
Britzen, S., Witzel, A., Krichbaum, T. P, et al. 2000, A&A, 360, 65 
Capetti, A. & Celotti, A. 1999, MNRAS, 304, 434 
Cohen, M. H., Lister, M. L., Roman, D. C, et al. 2007, ApJ, 658, 232 
Condon, J. J. & Mitchell, K. J. 1984, ApJ, 276, 472 
Conway, J. E. & Murphy, D. W. 1993, ApJ, 411, 89 
de Vaucouleurs, G. 1991, Science, 254, 1667 
Fan, J. H., Liu, Y., Yuan, Y. H., et al. 2007, A&A, 462, 547 
Fanaroff, B. L. & Riley, J. M. 1974, MNRAS, 167, 31P 

Gabuzda, D. C. 1995, Proceedings of the National Academy of Science, 92, 

11393 

Gabuzda, D. C, Pushkarev, A. B., & Cawthome, T. V. 2000, MNRAS, 319, 1109 
Ghisellini, G. 1993, Advances in Space Research, 13, 587 
Ghisellini, G., Padovani, P, Celotti, A., & Maraschi, L. 1993, ApJ, 407, 65 
Giovannini, G., Cotton, W. D., Feretti, L., Lara, L., & Venturi, T. 2001, ApJ, 552, 
508 

Gomez, J. L., Guirado, J. C, Agudo, I., et al. 2001, MNRAS, 328, 873 
Gomez, J. L., Marti, J. M. A., Marscher, A. P., Ibanez, J. M. A., & Marcaide, 

LM. 1995, ApL 449, L19+ 
Gong, B. 2008, MNRAS, 389, 315 

Gong, B., Kong, S. W., Xue, F, Li, Y, & Huang, Y. F 2011, ArXiv e-prints 
Hasinger, G. 1998, Astronomische Nachrichten, 319, 37 
Heidt, J. & Wagner, S. J. 1996, A&A, 305, 42 

Hewitt, A. & Burbidge, G. 1989, in A new optical catalog of QSO (1989), 0-+ 
Hook, 1. M., McMahon, R. G., Patnaik, A. R., et al. 1995, MNRAS, 273, L63+ 
Hopkins, P R, Hernquist, L., Cox, T. J., et al. 2006, ApJS, 163, 1 
Hough, D. H., Vermeulen, R. C, Readhead, A. C. S., et al. 2002, AJ, 123, 1258 
Hovatta, T., Valtaoja, E., Tomikoski, M., & Lahteenmaki, A. 2009, A&A, 494, 
527 

Jiang, L., Fan, X., Ivezic, Z., et al. 2007, ApJ, 656, 680 
Jorstad, S. G., Marscher, A. P, Mattox, J. R., et al. 2001, ApJ, 556, 738 
Karouzos, M., Britzen, S., Eckart, A., Witzel, A., & Zensus, A. 2010, ArXiy 
e-prints 

Kellermann, K. I., Lister, M. L., Homan, D. C, et al. 2004, ApJ, 609, 539 
Kellermann, K. 1., Sramek, R., Schmidt, M., Shaffer, D. B., & Green, R. 1989, 
AL98, 1195 

Krichbaum, T. P., Witzel, A., Standke, K. J., et al. 1994, in Compact 

Extragalactic Radio Sources, ed. J. A. Zensus & K. 1. Kellermann, 39 — h 
Kudryavtseva, N., Britzen, S., Witzel, A., et al. 2010, ArXiv e-prints 
Lahteenmaki, A. & Valtaoja, E. 1999, ApJ, 521, 493 

Laing, R. A. 1996, in lAU Symposium, Vol. 175, Extragalactic Radio Sources, 

ed. R. D. Ekers, C. Fanti, & L. Padrielli, 147-+ 
Lister, M. L., Cohen, M. H., Homan, D. C, et al. 2009, AJ, 138, 1874 



Lister, M. L. & Homan, D. C. 2005, AJ, 130, 1389 

Lister, M. L., Kellermann, K. I., Vermeulen, R. C, et al. 2003, ApJ, 584, 135 
Lobanov, A. P., Krichbaum, T. P., Witzel, A., & Zensus, J. A. 2006, PASJ, 58, 

253 

Lobanov, A. P & Zensus, J. A. 1999, ApJ, 521, 509 

Lowe, S. R., Gawrohski, M. P, Wilkinson, P N., et al. 2007, A&A, 474, 1093 

Madau, P, Ferguson, H. C, Dickinson, M. E., et al. 1996, MNRAS, 283, 1388 

Marscher, A. P & Gear, W. K. 1985, ApJ, 298, 1 14 

Mimica, P, Aloy, M., Agudo, I., et al. 2009, ApJ, 696, 1142 

Morganti, R., Oosterloo, T. A., Fosbury, R. A. E., & Tadhunter, C. N. 1995, 

MNRAS, 274, 393 
Nesci, R., Massaro, E., Maesano, M., et al. 2002, AJ, 124, 53 
Padovani, P 1993, MNRAS, 263, 461 
Padovani, P & Urry, C. M. 1991, ApL 368, 373 

Pearson, T. J., Browne, 1. W. A., Henstock, D. R., et al. 1998, in Astronomical 
Society of the Pacific Conference Series, Vol. 144, lAU CoUoq. 164: Radio 
Emission from Galactic and Extragalactic Compact Sources, ed. J. A. Zensus, 
G. B. Taylor, & J. M. Wrobel, 17-+ 

Pearson, T. J. & Readhead, A. C. S. 1981, ApL 248, 61 

Pelletier, G. & Roland, J. 1989, A&A, 224, 24 

Perucho, M., Lobanov, A. P., & Kovalev, Y. Y. 2009, in Astronomical Society of 
the Pacific Conference Series, Vol. 402, Astronomical Society of the Pacific 
Conference Series, ed. Y. Hagiwara, E. Fomalont, M. Tsuboi, & M. Yasuhrro, 

349-+ 

Perucho, M., Marti, J. M., Cela, J. M., et al. 2010, A&A, 519, A41+ 

Pollack, L. K., Taylor, G. B., & Zavala, R. T. 2003, ApJ, 589, 733 

Pushkarev, A. B., Kovalev, Y. Y, Lister, M. L., & Savolainen, T. 2009, A&A, 

507, L33 
Rees, M. J. 1966, Nature. 21 1, 468 

Roland, J., Britzen, S., Kudryavtseva, N. A., Witzel, A., & Karouzos, M. 2008, 

A&A, 483, 125 
Sol, H., Pelletier, G., & Asseo, E. 1989, MNRAS, 237, 411 
Spergel, D. N., Verde, L., Peiris, H. V, et al. 2003, ApJS, 148, 175 
Steffen, W., Zensus, J. A., Krichbaum, T. P, Witzel, A., & Qian, S. J. 1995, 

A&A, 302, 335 

Taylor, G. B., Vermeulen, R. C, Readhead, A. C. S., et al. 1996, ApJS, 107, 37 
Urry, C. M. & Padovani, P 1995, PASP 107, 803 

Vermeulen, R. C. 1995, Proceedings of the National Academy of Science, 92, 
11385 

Vermeulen, R. C, Britzen, S., Taylor, G. B., et al. 2003, in Astronomical Society 
of the Pacific Conference Series, Vol. 300, Radio Astronomy at the Fringe, 
ed. J. A. Zensus, M. H. Cohen, & E. Ros, 43-+ 

Vermeulen, R. C. & Cohen, M. H. 1994, ApJ, 430, 467 

Wagner, S. J. & Witzel, A. 1995, ARA&A, 33, 163 

Whitney, A. R., Shapiro, 1. 1., Rogers, A. E. E., et al. 1971, Science, 173, 225 
Witzel, A. 1987, in Superluminal Radio Sources, ed. J. A. Zensus & 

T. J. Pearson, 83-93 
Xu, Y.-D., Cao, X., & Wu, Q. 2009, ApJ, 694, L107 



